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Abstract 

The bi-Hamiltonian structure of the Benney hierarchy is revisited. 
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^^ . We show that the compatibility condition of the Poisson brackets pro- 

^^ , vides the genus zero free energy of a topological field theory coupled to 

2d gravity. We calculate the correlation functions via the Landau- 
Ginzburg formulation and derive the string equation based on the 
twistor construction. Moreover, by using the approach of Dubrovin 
and Zhang, we compute the genus one correction of the Poisson brack- 
ets and compare them with the Oevel-Strampp's brackets of the Kaup- 



_C^ \ Broer hierarchy. 



I. INTRODUCTION 



In the past decade, the developments of integrable systems have made many 
important influences on theoretical physics and pure mathematics. Among others, 
those works concerning the relationship to topological field theories (TFT) and string 
theories have been paid much attention in the frontier subjects (see for a review). 
In particular, Witten |^0|] and Kontsevich |^ show that the partition function of 
2d topological gravity is equivalent to a particular tau-function of the Korteweg-de 
Vries hierarchy characterized by the string equation. Now it is generally believed 
that 2d TFT coupled to 2d gravity can be formulated as integrable hierarchy of 
nonlinear partial differential equations. 

In general, 2d TFT can be classified by the solutions of the Witten-Dijkgraaf- 
E.Verlinde-H.Verlinde (WDVV) equations of associativity [^,§] in the sense that a 
particular solution of WDVV equations provides the primary free energy of some 
topological model. In fact, various classes of solutions to the WDVV equations 
have been obtained (see ||31| , p!0|JTT[| and references therein), which turn out to be the 
tau-functions of dispersionless integrable hierarchies. Accordingly, investigating the 
solution space of the WDVV equations will deepen our understanding of 2d TFT. 

It is well-known that Poisson structures of dispersionless integrable hierarchies 
have the form of hydrodymanic type |M. Due to this fact, the integrability (bi- 



Hamiltonian structure) of the WDVV equations can be formulated in a geometric 
way called Frobenius manifolds ]ll|. Based on this geometrical construction, the 
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higher genus extension of the Poisson structures [0 and Virasoro constraints 
of the associated integrable hierarchies are given. 

The purpose of this paper is to study a dispersionless nonstandard Lax hierarchy 
from the TFT point of view, which is a modification of the dispersionless Kadomtsev- 
Petviashvili (dKP) hierarchy p8| , |38 
form 



The Lax operator we would consider has the 
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which satisfies the hierarchy flows (Ti 
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{^>f,^}, 



1.1^ 



;i.2) 



where L^i means the polynomial p" + ■■■ + (■ ■ ■)p, i.e., we cut off the terms after p 
of the expansion (in p) of L" and the Poisson bracket { , } is defined by |0 

df dg df dg 



{f{p,X),g{p,X)} 



;i.3) 



dp dX dX dp 

We remark that the Lax operator ( |1 . 1| ) can be obtained from that of the dispersion- 
less modified KP (dmKP) hierarchy p4|,3 via truncations. Thus the nonstandard 
Lax hierarchy ( |1.2|) is referred to the constrained dmKP hierarchy [||. 



In this paper, for simplicity, we shall concentrate on the Lax operator of the form 

L = p + v^ + v'^p~'^ (1.4) 



which satisfies the nonstandard Lax equations 

dL 



dTn 



{Ll,.L} 



1.5) 



(we will leave the cases for A^ > 2 in a subsequent publication). The first few flows 
are 
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(v^r + 2v^ 
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where the simplest equation (r2-fiow) being the Benney equation describes long 
waves in nonlinear phenomena (here we only consider the 1 + 1 dimensional reduction 
of the Benney system in 0). The whole equations ( |1.5|) form what we call the 
Benney hierarchy which has been intensively studied during the past two decades 
(see, for example, |p2| , ^|25| , ^|33| , p| and references therein). 

Remark 1. The algebraic and Hamiltonian structure associated with the kind of 
Lax operators 



1.1 



have also been investigated in [^ where the Lax equations are 
defined by the bracket {A, B} = pdA/dpdB/dx — pdA/dxdB /dp with respect to 
the decomposition A = A>o © A<o for the pseudo-differential operator A and the 
relationship to the dispersionless Toda hierarchies (dToda) [^,3S| is established as 
well. 

Recently, the bi-Hamiltonian structures associated with the Benney hierarchy 
( |1.5[ ) have been investigated based on the theory of classical r-matrices |]3^. From 
that one can associate with a free energy coming from a 2d TFT. In particular, by 
this free energy, we construct an additional hierarchy generated by Hamiltonians 
with logarithmic type which together with the ordinary hierarchy are identified as 
flows in genus zero TFT coupled to 2d topological gravity. The basic idea is that 
the dispersionless Lax operator (|1.4]) can be viewed as a superpotential in Landau- 
Ginzburg (LG) formulation of TFT [p9i 3- Thus, according to the LG theory, the 
variables v^ and f ^ are identified as the fundamental correlation functions and their 
dynamic flows turn out to be the genus zero topological recursion relations p0|,|9[| of 
the associated TFT. 

Moreover, in order to establish the string equation describing the gravitational 
effect, we construct the twistor data for the Benney hierarchy by using the Orlov 
operator corresponding to the dmKP hierarchy Q. We show that a remarkable 
feature of Benney hierarchy is the flows generated by the additional logarithmic 



Hamiltonians can be expressed by the logarithm of the Lax operator and are well- 
defined only after suitable constraint included in twistor data. 

Finally, from the genus zero free energy of the Benney hierarchy, we compute the 
associated G-function to construct the genus one free energy and then quantize the 
Poisson brackets of the Benney hierarchy using the Dubrovin-Zhang's (DZ) approach 
to bi-Hamiltonian structure in 2d TFT |l^. On the other hand, we also "quantize" 



these Poisson structures from the Oevel-Strampp's (OS) brackets of the Kaup-Broer 
(KB) hierarchy p9| , |35[| . We find that after appropriate differential substitutions, they 



are matched up to genus one correction. 

The paper is organized as follows. In next section we compute the primary free 
energy from the bi-Hamiltonian structure of the Benney hierarchy and then intro- 
duce the additional logarithmic fiows commuting with the ordinary Benney fiows. In 
Sec. Ill we compute topological correlation functions using the LG formulation and 
show that the Benney hierarchy can be derived from genus zero topological recur- 
sion relations. Sec. IV is devoted to finding the twistor data to establish the string 
equation including the additional logarithmic fiows. In Sec. V we show that the 
genus one correction of the Poisson brackets obtained by DZ quantization coincide 
with the OS brackets after appropriate differentiable substitutions of the dynamical 
variables. In the last section we discuss some problems to be investigated. 

II. BI-HAMILTONIAN STRUCTURE AND FREE ENERGY 

In this section, we shall investigate the relations between the bi-Hamiltonian 
structure and its associated free energy of the Benney hierarchy. It turns out that 
we can find additional logrithmic hierarchy using the free energy, which together 
with the ordinary hierarchy ( |1 . 5|) form the fiows of gravitational couplings in TFT 
of genus zero, the topic to be discussed in next section. 

The bi-hamiltonian structure of the Benney hierarchy ( |1.5| ) is given by @,^ 

dv 6Hn+i _ 6Hn 
•Ji — ; — o'2" 




/ 2d dv^ 



la „,2A I a„,2 I ■ (2-1; 



V d V d + dv 



with Hamiltonians defined by 



Hn = - resL'' 
n J 



where d = d/dX and resL" being the coefficient of p ^ of L". We list some of them 
as follows: 



Hi 
H2 
H3 
H, 



v\\ 






Remarks 2. The second bracket J 2 in fact reveals the classical limit of Virasoro- 
t/(l)-Kac-Moody algebra P| with v^ being the Diff5'^ tensor of weight 2 and v^ 
a tensor of weight 1. This is due to the fact that the DiffS*^ flows are just the 
Hamiltonian flows generated by the Hamiltonian li\ = Jv^. 

Besides the concept of integrability, the geometrical means of the Poisson brack- 
ets (|2.1| ) is profound. The essential idea is based on the fact that the bi-Hamiltonian 
structure Ji and J2 can be written as 



Ji 
J2 




g°^^[v 



)9 + rf( 



V]V 



X- 



where F"^! 



v) is the contravariant Levi-Civita connection of the contravariant flat 
metric g°'^{v). Both Ji and J2 are Poisson brackets of hydrodynamic type introduced 



by Dubrovin and Novikov |T^ . The existence of a bi-Hamiltonian structure means 
that Ji and J2 have to be compatible, i.e, J = Ji + XJ2 must be a Hamiltonian 
structure as well for all value of A. The geometric setting of this bi-Hamiltonian 
structure of hydrodynamic system is provided by Frobenius manifolds [|TT],|12[. One 
way to define such manifolds is to construct a function F{t^,t'^, ■ ■ ■ ,t"^) such that 
the associated functions, 



^a/37 



d^Fjt) 



(2.2) 



satisfy the following conditions ||11|| : 

• The matrix rjap = ciap is constant and non-degenerate, (for the discussion of 



degenerate cases, see [BH]) 



The functions c^^ 



define an associative commutative algebra with a 
unity element. The associativity will give a system of non-linear PDE for F{t) 



V°"c^f3-/ 



d^F{t) 



:v 



A/1 



d^F{t) d^F{t) 



dt^dt^dt^ ' dtf'dPdt'' dt^dt^dt^ ' dtf^dt^dt^ 



\fi 



d^F{t) 



(2.3) 



The functions F satisfies a quasi-homogeneity condition, which may be ex- 
pressed as 

CeF = dpF + (quadratic terms), 
where E is known as the Euler vector field. 



Equations (|2.3|) constitute the WDVV equations [pjipj arising from TFT (see Sec. 
III). A solution of the WDVV equations will be called primary free energy. Given 
any solution of the WDVV equation, one can construct a Frobenius manifold M. 
associated with it. On such a manifold one may interpret r/"^ as a flat metric and 
t" the fiat coordinates. The associativity can be used to defines a Frobenius algebra 
on each tangent space T^M.. This multiplication will be denoted hy u-v. Then one 
may introduce a second flat metric on M. defined by 

g""^ = E{dr-dt^), (2.4) 

where dt" ■ dt^ = d^^dV = rj'^^c^ dV . This metric, together with the original metric 
r;"^, define a flat pencil (i.e, 77"^ + Af?"^ is fiat for any value of A). Thus, one 
automatically obtains a bi-Hamiltonian structure from a Frobenius manifold M.. 



The corresponding Hamiltonian densities are defined recursively by the formula |TT 



where n > l,a = 1,2,- ■ ■ ,m, and h'^^ = rfajst^ ■ The integrability conditions for this 
systems are automatically satisfied when the c^^ are defined as above. 

To find the free energy associated with the Benney hierarchy ( p..5|) , one might 
set t^ =v^ = h2 , t^ = v"^ = hi and thus 

r/nt)=f?il, ^"^(t)=f'tV (2-6) 



We remark that the fiat metrics ri°'^{t) and g°'^{t) satisfy ri°'^{t) = dg"^{t)/dt^ and 
J t^ and /t^ turn out to be the Casimirs for the bi-Hamiltonian structure of the 
hierarchy. In fact, those c^^ can be determined by ( p.5|) and ( p.6|) . For a = 1 we 
have 

M-) = ^^ = l/resL", ihf^=t-) 

[n — 1)! nl J 

which, up to a normalization, are the Hamiltonian densities of the Benney hierarchy 
and 

•^ll ~ '^22 ~ U, C21 = C12 = 1- l^- ' J 



Then, from ( p.TD and (|2.2| ), we get immediately the free energy 



1 j.2\ -'-/j1\2j_2 I '- u2\2 fi ,j.2 "^ 



F(t\t^) = -(m^ + -(t^r(^logt^-2 



(2.8) 



Also, from {UM, 



and 



we can obtain the associated Euler vector field 



E = t 



1 d 



+ 2r 



d 



which implies the quasi-homogeneity condition: 

CEF{t) = AFit) + iff. 



Next, we turn to the hierarchy corresponding to a = 2 with /i2 = i^ = v^ . 



M 



Using dU) and (|j), we get 






(tT 

6 



+ t^(logt^-l), 
+ tV(logt'-l) 



h^' = ^ + \{t'ft\\oge - 1) + \{ef (logt^ - ^) , 



/zi^) 



24 
(t^^ 1 
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+ i{t'fe{\oge-i) + \t\ef{\oge 



Motivated by the work of ||2T|,[T8|, the Hamiltonian densities /i2 can be expressed 
as 

2 



h^^ 



with the prescription 

logL = log(p + t^ + t V^) 



n\ 



-res[L''(logL-c„)], 



t' 



= - log t' + - log(l + t^p-^ + ep ') + - log ( 1 + -p + -p" 
and c„ = I]?=i ^, Co = 0. Then the Lax flows corresponding to /i2 are 



(2.9) 



-j=i i 






2{5„,L}, 5„=[L"(logL-Cn)]>i 



(2.10) 



or, in terms of bi-Hamiltonian structure 



Jl ? — J2 



dTn 



5v 



(5v 



where the Hamiltonians Hn are defined by 

^„ = '^jm^[L^{\ogL-c^)] = {n-l)\Jh 



These Hamiltonians generate additional fiows which are compatible with the ordi- 
nary Benney fiows. We will see later that it is these logrithmic fiows ( ^.101 ) which 



together with the ordinary fiow ( |1.5|) implies that the Benney hierarchy can be 



formulated as a 2d TFT coupled to gravity. 



III. TOPOLOGICAL STRING AT GENUS ZERO 

In this section we would like to set up the correspondence between the Benney 
hierarchy and its associated TFT at genus zero. Let us first recall some basic notions 
in TFT. 

A topological matter theory can be characterized by a set of BRST invariant 
observables {Ci, O2, •••} with couplings {T"} where Ci denotes the identity opera- 
tor. If the number of observables is finite the theory is called topological minimal 
model and the observables are refereed to the primary fields. When the theory 
couples to gravity, a set of new observables emerge as gravitational descendants 
{an{Oa),n = 1,2, ■ ■ ■} with new coupling constants {T"'"}. The identity operator 
Oi now becomes the puncture operator P. For convenience we can identify the 
primary fields O^ and the coupling constants T" to ao{Oa) and T"'°, respectively. 
As usual, we shall call the space spanned by {T'^'^-jU = 0, 1,2, ■ ■ ■} the full phase 
space and the subspace parametrized by {T"} the small phase space. These cou- 
pling times decribe the perturbative fiows with respect to the corresponding critical 
theory (in which T""'" = 0). 

For a topological model the most important quantities are correlation functions 
which describe the topological properties of the manifold where the model lives. The 
generating function of correlation functions is the full free energy defined by 

00 00 

HT) = E-^.(^) = E(e^-"^'^'""^''"^). (3-1) 

9=0 g=0 

where (■ ■ ■)g denotes the expectation value on a Riemann surface of genus g with 
respect to a classical action. In the subsequent sections, we will omit the exponen- 
tial factor without causing any confusion. Therefore a generic ?7i-point correlation 
function can be calculated as follows 

(cr„^(0„Jcr„2(C„2) ■ ••(T„„(0„„))g = ^y„j „^^y„2^„2 _ _ _ 57-a™,n„, • '^3-2) 

In the following, we shall restrict ourselves to the trivial topology, i.e., the genus zero 
sector {g = 0) since this part is more relevant to dispersionless integrable hierarchies. 



In particular, the genus zero free energy restricting on the small phase space is the 
primary free energy defined by 



T, 



T'^^i'^ 2^a,n>l^Q 



F{t). 



(3.3) 



Let us define some genus zero correlation functions on the small phase space. 
The metric on the space of primary fields is defined by 



{PO^Oi,) = f]^^ 



(3.4) 



When r]ai3 is independent of the couplings we call it the flat metric and the couplings 
T" the fiat coordinates. In fact, a three-point function in the small phase space can 
be expressed as 



(Oo^OfsOj) 



OF 



Caf3^ 



(3.5) 



dT^dTf^dTy 
which provide the structure constants of the commutative associative algebra 

Oo^Op = clp, (3.6) 

with constraints c^^ = rp^Caap, ciap = rjap. The associativity of c^^^, i.e.. 






will give the WDVV equations (|2l| ). 

Now, Let's return to the Benney hierarchy. Since the Benney hierarchy is a 
two- variable theory, thus only two primary fields {Oi = Op = P, O2 = Oq = 
Q} are involved in the TFT formulation and we shall identify f^|T<».">i=o = T^ 
and v'^\xo,n>i^Q = T^ on the small phase space. Therefore the Lax operator in 
small phase space is written as L[z) = z -\- T^ -|- T'^z~^ which can be viewed as a 



superpotential in LG formulation of TFT ||3^J§]. According to the LG theory, the 
primary fields are defined by 



Op(z) 



dL{z) 



Ooiz) 



dL{z) 



~-i 



(3.7) 



QJ^P ' ^vv y QJ.Q 

which can be used to compute the three-point correlation functions through the 



formula ||39|,|8[|: 



Ca/s-f — resL/=o 



Oa{z)Op{z)0,{z) 
dMz) 



(3. 



It is easy to show that ( |3.8| ) reproduces the previous Cafs^ and F on the small phase 
space. In particular, the fiat metric on the space of primary fields is given by 



Vpq = Vqp 



Vpp = Vqq = ^ 



as we obtained previously. Now we can impose the fundamental correlation functions 
as {PP) = d^F/dirPy = T«, {PQ) = d^F/dT^dT^ = T^, and {QQ) = logT«. 
Although these two-point correlation functions are defined on the small phase space, 
however, it has been shown [^ that they can be defined on the full phase space 
through the variables v^ and f ^ in which the gravitational couplings T"'" do not 
vanish. Hence it is easy to write down these genus zero two-point functions on the 
full phase space and obtain the following constitutive relations: 



{PP) 

(PQ) 
{QQ) 



d{Tpy 

dTPdTQ 
log(PP) 



(3.9) 



which will be important to provide the connection between the Benney hierarchy 
and its associated TFT. 

Remark 3. In general, for the two-primary models, {QQ) = f{{PP)) where the 
function f{x) is model dependent 0. In 0, the same relation f{x) = logx has 
been imposed to extract the nonlinear Schrodinger hierarchy from the hermitian 
one-matrix model at finite N. However, for the CP^ model |^,|9[|, f{x) = e^. 

Based on the constitutive relations ( p.9|) , we can identify the gravitational flows 
for v^ and v^ in the full phase space as the Lax flows by taking into account the 
genus zero topological recursion relation |E0l3: 



{ar,{0^)AB)= Y. n{a^^^{0^)Op)r^f'^0,AB), 
For example, setting ra = 1, Oa = P and A = P, B = Q then 

{<yi{P)PQ) 

{PP){QPQ) + {PQ){PQQ) 
{PP){QQ)'+{PQ){PQ)' 



a 



P,Q. 



(3.10) 



dTP' 



where we denote /' = df/dT^ = df/dX . Similarly, taking A = P, B = P we have 



On the other hand, taking ra = 2 we get 

dv^ ri 



{ai{P)PP) = {v^v 



1„.2n/ 






1„,2 



{v^y + 2v^v 



2n2 



(^i)^^2 ^ (^.) 



(3.11) 



Comparing the above equations with the Lax flows (|1.6| ), we shall identify T„ = 
TP'^-i/n^ (n = 1,2,---). 

Next, let us turn to the T*^'" flows. Choosing Oa = Q and using the topological 
recursion relation ( 3.10 ), we obtain 






i(v'r + v\\ogv' - 1) 

(t;^)2logi;2 + 2t;2(logw2_2^ 



,1n3 



+ 2t;^t;^(logt;^-l) 



(3.12) 



which are nothing but the Lax flows associated with the logarithmic operator i?„ 
for n = 1,2 under the identification T„ = T*^'". 

It turns out that equations ( p.ll| ) and (|3.12| ) can be recasted into the Lax form 
in terms of coupling times T°'" as 



dL 
dL 



{Ll\\L}, 



n + 1 
2{[L"(logL-c„)]>i,L}. 



(3.13) 



and the associated commuting Hamiltonian flows with respect to the bi-Hamiltonian 
structures are thus given by 



9v 

Q^P,n 



{-f?^P,n+l, v}i = {ifp^rt, v}2, Hp^n 

{-f^Q,n+l, V}l = {iJQ,„, V}2, Hq^ri 



lesL' 



n+l 



n{n + 1) 
^y'res[L"(logL-cO] (3.14) 



n 



where n = 0, 1, 2, ■ ■ •. 



Remark 4. We notice that ( |2.8|) corresponds to the primary free energy of the dKP 
system in the iyo,i-model |^ defined by the Lax operator L = p + v^i/{p — s) under 
the identification s = T^,v^i = T^. They also show that another choice leads to 
the dToda system in the CP^ topological sigma model |pll , [l8 |. 



Furthermore, using the constitutive relations ( p.9| ) and the Lax flows (|3.13| ), we 
have 



dv' d{an{0^)Q) 



grpa,n grpP 



Q^a,n 



QrpP 



\-'~^a,n) ) 



where R'^}^ are the analogues of the Gelfand-Dickey potentials [^] of the KP hierar- 
chy, given by 

R% = MP)Q) = ;^(^"^')o, 

RqI = MQ)Q) = 2[L"(logL - c„)]o, 
R^Z = MP)P) = ^^resL'^+S 

RqI = K{Q)P) = 2res[L"(logL - c„)]. (3.15) 

For example, the two-point correlators involving the first and the second descendants 



are 



{a^iP)P) = vh', 
{a,{P)Q) = ^ivr + v\ 

{a^{Q)Q) = vHogv^ 
{a,{P)P) = iv'fv' + (vY, 

{a2{Q)P) = l{v'f + 2v'v\\ogv'-l), 
(^2(Q)Q) = (v^yiogv^ + 2t;2(logt;2 _ 2). 

Finally, we would like to remark that the last two equations of (|3.15|) can be 



integrated to yield 



^ "(2) 1 .o.r™+2 



(a„(P)) = ^R^+i = r^TTT^^"^'^ 
n + 1 {n + lj[n + 2) 

MQ)) = -^RqI+i = ^res[L"+i(logL - c„+0]- 

which are just the LG representation for one-point functions of gravitational decen- 
dants at genus zero. 



IV. THE TWISTOR DATA AND STRING EQUATION 

In this section we would like to discuss the string equation of the 2d TFT as- 
sociated with the Benney hierarchy, which govern the dynamics of the variables f " 
(or fundamental correlators) in the full phase space. The Lax formulation in Sec. 



II for Benney hierarchy is in fact similar to the formulation of the dToda type hi- 
erarchy [^,^. Based on this observation, we can reproduce the Benney equations 
by imposing constraints on the Lax operators and the associated Orlov operators 
of the dmKP hierarchy 0] through the twistor data (see below). We shall remark, 
however, that the flow equations corresponding to /i2 are absent in the standard 
formulation of the dToda type hierarchy. So we have to properly extend the stan- 
dard Orlov operator to include the additional hierarchy equations ( p.lO| ). We will 
follow closely that of p^ to show that the constraints imposing on the twistor data 



implies the string equation of the Benney hierarchy. 

Let us consider two Lax operators fi and fi with the following Laurent expansions 

(Ti = xy. 

oo 

;, = p + ^^;-(r,rK", 

n=0 

oo 

/i-1 = iJo{T,f)p~' + E Sn+i(T,f K- (4.1) 

n=0 



which 


satisfy the 


commuting 


Lax flows 




dTn 


= {Bn 


,M, 


dfn 




dfi 


- \U 


n\ 


dfi 



{i?„,/i}, 

.^ .^.,.., .^ {Bn,-^}, (ri = 1,2,3 ■••) (4.2) 

where the Poisson bracket {, } is defined as before and 

Bn = (/^")>i, Bn = (/i"")<o. 

Next, we consider the Orlov operators corresponding to dmKP 

oo oo 

M = E riT„/i"-i + E ^n(T, f )/i-" 

ra=l ra=l 

oo oo 

M = - E nf^r""'^ + X+Y. MT, f)fi^. (4.3) 

n=l n=l 

with constraint 

{/i,M} = l, {fi,M} = l. (4.4) 

In fact, the coefficient functions Wn and Wn in the Orlov operators are defined by 
the above canonical relations and the following flow equations 

dM dM 

^ = {Bn,M}, -^ = {B^,M}, (4.5) 

dM ~ dM ~ ~ 

-— = {5„,M}, ={i?„,M}. (n= 1,2,3 ■■■) (4.6) 

CJ n dT„, 



Inspired by the twistor construction (or Riemann-Hilbert problem) for the solution 
structure of the dToda hierarchy ||37| . |38[| , we now give the twistor construction for 
the Benney hierarchy. 

Theorem 1 Let f{p, X),g{p, X), f{p, X),g{p, X) be functions satisfying 

{f{p,X),g{p,X)} = 1, {f{p,X),~9{p,X)} = 1. (4.7) 

Then the functional equations 

/(/x, M) = /(/2, M), gifx, M) = ~gifi, M) (4.8) 

can get a solution of U-W o'^^ UA )- ^^ call the pairs {f,g) and {f,g) the twistor 
data of the solution. 

The proof is provided in appendix A. 

To reduce the above Theorem to the Benney hierarchy, we have to impose the 
following constraint on the Lax operators 

L = /i = /i-\ (4.9) 

that is, f{p,X) = p and f{p,X) = p~^. As a result, the time variables T„ can be 
eliminated via the following identification: 

T = -T 

From ( [4.7| ), the twistor data g{p, X) and g{p, X) can be assumed the following form 

oo 

gip,X) = X - ^ nT„p"-\ ~gip,X) = -Xp\ (4.10) 

n=2 

where the second part of g{p, X) is responsible for the string equations (see below). 
By the Theorem 1 and equation (|4.1CI|) , we get the following constraint for the Orlov 
operators: 

oo 

M-Y, nTnfi'''^ = -fi^M. (4.11) 

n=2 

It's the above constraint that leads to the string equations. 

So far, the twistor construction only involves the Lax flows ( |1.5| ). To obtain the 
string equations associated with the TFT described in Sec. II, we have to modify 
the Orlov operator M to include the additional flows ( [2.101 ). Namely, it's necessary 
to introduce the flows generated by the logarithmic operator 

5„ = [L"(logL-c„)]>i. 



where we have imposed the constraint (|4.9|) and used the prescription of the series 
expansion ( p. 91) for logL. Let T„ be the time variables of additional flows generated 
by Bn then the Orlov operator M is deformed by these new flows to M' so that 



dM' 



2{Bn,M'}. 



(4.12) 



To construct the modified Orlov operator M' , it is convenient to using the dressing 
method |^8|]. Let us first express the original Lax operator /i and its conjugate Orlov 
operator M in dressing form (similarity transformation) 



/i = e^d®(p), 



M 



^ade 



Y^ nTnP 



n-l 



(4.13) 



Kn=l 



where 



e^<^f(g)=g + {f,g} + -{f,{f,g}} + .... 

One can understand that this is the canonical transformation generated by 0(T, T, p) 
(for the T-dependence, see below) and its flow equations (Sato equations) can be 
written as [BSIl 



VT„,ee = i?n-ead®(p" 



(4.14) 



where 



^-^^?„(^ 



(adO)^ 






It is easy to show that (|4.14| ) together with (|4.13|) implies the flow equations ( ^75|) . 
In contrast with equation ( 4.14| ), the T„ flows for are given by 



Vf„,e0 = 25„ - e^d®[p"(logp - c„)]. 
and a similar argument reaches the modified Orlov operator 



M' 



,ade 



£ nT^p"-^ + 2 £ nf„p"-i(logp - c„_i) 



Vra=l 



n=l 



^ nT,,lx^-' + Y, w„/i-" + 2 Y ^T„/i"-i(log/i - c„_i) 



n=l 



n=l 



n=\ 



(4.15) 



which satisfies the additional flow equations ( |4.12| ). 

Now we are in the position to derive the string equation. Following |2^ we 
decompose M' into the positive power part (M')>i and the non-positive power part 
(M')<o by using equations ( |4.11| ) and (|4.15|) . It turns out that 



M'-^nT„/i"-M =2 

V n=2 / >i 



^n— 1 J 



71=1 



.n-1 



V n=2 / <o 



^nf„/i"-^(log/i-c, 

n=l 

2 ^ nfnBn-i 
=1 

(/i^M) 

oo 



'<0 



-I <0 



\n=l 



<0 



where, by the definition ( |4.1| ), we have used the fact that (/i ")>i = (/i")<o = for 
n> 1. Hence 



.n-l 



n=l \n=l / <o 



(4.16) 



which together with the canonical commutation relation {g{L^ M'), L} = —1 implies 

-1. (4.17) 



^ dL ^ - dL 

n=2 O-l-n-l n=l C^n-l 



After shifting Ti — ^ Ti — 1 and making the identification T„ = T^'" ^/n and T„ 
T*^'*^ as described before, we have 



/-^ ffpP.n-l /-^ 

n=l ^-^ n=l 



9T^ 



gXQ,n-i ■ 



(4.18) 



Now taking the zero-th order term of the above equation and using the constitutive 
relation ( p. 91) yield 



^^ = i + E E nT-"(a„_i(o„)pg). 



(4.19) 



n=l a=P,Q 

By integrating the above equation we obtain the universal string equation at genus 



zero 



t;^(T) = T^+E E riT"'"(a„_i(a)Q), 

n=l a=P,Q 



V 



■(T)=T« + E E riT°'"(a„_i(0«)P) 



(4.20) 



n=l o=P,Q 

which describe the dynamics of v"' associated with the gravitational background. 



V. GENUS ONE CORRECTION OF POISSON BRACKETS 



In this section, we will compute the genus one correction to the bi-Hamiltonian 
structure of the Benney hierarchy by using the DZ approach to integrable hierarchy 
associated with TFT |]TB[ , which consists of the following two main ingredients: 



introducing slow spatial and time variables scaling 



rpc 



er 



n = 0,1,2, 



(5.1^ 



changing the full free energy as 



fl=0 



9' 



where e is the parameter of genus expansion. Thus all of the corrections become 
series in e. 



To get a unambiguous genus one correction of the Hamiltonian flows ( |3.14]) one 
may expand the flat coordinates up to the e^ order as 



.2^(1) 



ta = tr+eX'^ + 0(e^), t^ = Vapt^ 



(5.2) 



where t^^^ is the ordinary Benney variables Va satisfying (|3.14|) and t^^^ is the genus 
one correction defined by 



t 



(1) 



dT^dX 



(5.3) 



Then there exists a unique hierarchy flows of the form |15 



dTP,n 



{r(X),/J^,„+i}i + 0(e^ 
{r(X),i7^,„}2 + 0(e^) 



with 






(5.4) 



1,2 
(5.5) 



That means under such correction the Poisson brackets Ji and J2 and the Hamil- 
tonians will receive corrections up to e^ such that the Hamiltonian flows ( ^.4|) still 
commute with each other. 



In IjT^ , the genus-one part of the free energy has the form 



^i(T) 
where the matrix M^ is given by 



-logdetM^{t,dxt) + Git) 



t=v(T) 



(5.6) 



M|(t,5xt) =c%^{t)dxt\ 



and G{i) is a certain function satisfying the Getzler's equation p^. The first part of 
the formula (|5.6|) is quite simple on the small phase space, whereas, the second part 
describes, in the topological sigma-models, the genus one Gromov-Witten invariant 
of the target space and satisfies a complicated recursion relations [Q. For the 
primary free energy ( |2.(j| ), the G-function satisfies the following simple ordinary 
differential eqaution [p!5|] : 



48 (t 



2\2 



d^G 



o. 2dG 
2Af-— 



which can be easily solved as (up to a constant) 

1 



Git) 



12 



logt^ 



(5.7) 



The G-function of the Frobenius manifold also satisfies the quasi-homogeneity condi- 
tion CeG = —1/6. The above G- function can also be derived using the tau function 
of the isomonodromy deformation problem arising in the theory of WDVV equations 
of associativity [0. We briefly describe the derivation in appendix B. 
On the other hand, a simple computation yields 



M° 



^X ^X 



which together with 



implies 



^i = ^lo^ 



{t\fe - {t\f] -\\ogt\ 



Using dp^ and Ti and consulting the procedure developed by Dubrovin and Zhang 
(c.f. Theorem 1, Theorem 2 and Proposition 3 in [^]), after a straightforward but 
tedious computation, we obtain the genus one correction of the first Poisson bracket: 



{ti(X),ti(F)K = + 0(6^), 
{t\X),t\Y)}, = 5'{X-Y) + 



6 



-5"{X - Y) 



{t\X),t\Y)}, = 5\X-Y) + ^- 



0|2 

^*^ 5"{X 



{e{X),t\Y)], 



+ O(e^ 




^5"\X-Y) 



y) 



+ 0{e' 



+ 0(e^ 



(5.^ 



On the other hand, for the second bracket, we get 



{t\X),t\Y)}, = 25'{X-Y) + - 



it')' 
2 



/-2 \2' 



wj '^'' - '■* 



+ |^^"(^ -Y)- f.nx - Y) 



0(e' 



{t\X),t'iY)}2 = t\5{X -Y)+ t'5'iX - Y) 



e^ 
+¥ 



tH^x 



t 



e 



{t\X),t\Y)}2 = t'6'{X-Y) 



f U"{X-Y)--6"\X-Y) 



6 



t txx _|_ 2tx^x 



2t^{t 



X) 



(t 



2^2 



it 



2^2 



it 



2\3 



+ 



'2tH\ 



2t 



t 



aiV\ +2/VM —4-2 SI V VA , 0+2x'/\^ VA , rM A' 



{t\X),t\Y)}2 = tj,6{X -Y) + 2t'6'{X -Y) + O(e^). 



t' ) 



5'(X - Y) 



(5.9) 



Also, we can derive the genus-one corrections of the Hamiltonians, ii^^\ and some 
of them are: 






E 



(1) 

P,2 



o-{l) 



i/, 



(1) 



6t2 i 



(t 



1 \2 



t tj^tj^ 



6 

■-Xy 



12t2 ' 24(t 



+ 



6t2 



itY' 



2^2 



, (41! 



dX, 



t\t^'^ 



'^Xj 






rfX. 



Therefore the corresponding commuting flow equations up to genus-one corrections 
are given by 



t' 



t' 



Tpfi 



^Q,0 



Tq, 



t' 



0(t' 



X 



{t^f 



tH 



logt^ 



U2 



2*xx 



3{t|) 



2 \2 



(t2)2 



X 



24 I zi/i - 4fjytj, 



0(e^ 



X 



(t2)2 



e 

t^ ;^ + 24 

tMogt^ 

m' + tHiogt' - 1] 



2*2^ 2(t^ 2(t2j2 



(t2)2 





(t2)3 



X 



+ 0(6^), 



X 






+ 



t2 



24 



^^ ''XX 
(t2)2 

4+2 

'±1- ^^ ^ 



+ 



(t2)2 



+ 



(f)3 



5(4) ^ 

(t2)2 



(t2)2 



+ 0(e^ 



X 



Finally we would like to show that the genus one correction of the Poisson brack- 
ets ( |5.8| ) and ( |5.9| ) can be rederived from the "quantum" brackets associated with 
a dispersive counterpart of the Benney hierarchy. That means the loop correction 
can be viewed as the dispersive effect of the hydrodynamic type Poisson structure. 
Let us consider the "quantum" Lax operator of the form 

K = ed + u^ + {eOy^u'^, d = d/dX 



which is just the Lax operator of the KB hierarchy discussed in [^,0 under the 
scaling ( ^.1|) . The bi-Hamiltonian structure associated with K has been obtained 
by Oevel and Strampp p5| as follows 



{I, Jh = J 



res 



SK \5K. 



1,2 



where / and J are functionals of K and the Hamiltonian maps Qi are defined by 



Qi : 



n. 



51 

Ik 

6I_ 

Ik 



'5V 
6K 



>i 



K- 



+ 



5r 

r-X 



K-K 
SI 



51 



K 



>-i 



res 



5K 



.K 



5K^ 



,K 



-K\ - 



K 



5r 

5K 



.K 



Jk'- 



.K 



K 



with 



— - — av^— 



Using the Hamiltonian flows dK/dt = Qi{6H/SK), we can easily read off the "quan- 
tum" Poisson brackets 

{u\X),u\Y)}, = 0, 
{u\X),u\Y)}, = 6'{X~Y), 
{u'{X),u\Y)}, = 5'{X-Y), 
{u\X),u\Y)}, = 



for the first structure and 



{u\X)y{Y)}, = 25\X-Y), 

{u\X),u\Y)}2 = u],5iX -¥)+ u'5'iX -Y) + eS"iX - Y), 

{u\X),u\Y)}2 = u'6'{X -Y)- e6"{X - Y), 

{u\X), u\Y)}2 = uj,6{X -Y) + 2uH\X - Y) 

for the second. As a result, the first structure gets no correction, whereas the second 
structure receives a first order correction. So far, everything is exact. However if we 
define the following substitution for the fiat coordinates t": 



t\T) = u^ - eilnu" 



u 



IX 



X 



24 V m2 



+ 



X 



72 



{\nu^)xxu^ -{u]c 



1 ^2 



[u 



2^2 



O(e^), 



X 



V ^ 2 ^ 8 ^ ^ 144 I m2 



Oie' 



(5.10) 



XX 



where the right-hand side of t" is constructed from c}*f°/(9(T'^)*|„Q=„Q, then a 
straightforward but lengthy calculation shows that the first and the second Poisson 
brackets for t" coincide with equations (|5.8|) and ( ^.9|) modulo 0[t^\ Furthermore, 
using ( p.lO| ), it is not hard to check that the dispersion expansion of the Hamiltoni- 
ans and hierarchy flows defined by the KB hierarchy and those defined by Hpn and 
T^'"-flow coincide, modulo O(e^). In this sense, the parameter e of genus expansion 
characterizes the effect of dispersion. 
Remarks 5. The free energy and G- function associated with the bi-Hamiltonian 



structure of Drinfeld-Sokolov reduction of Lie algebra B2 are ||T^ 



;(^ 



U2^2 



15 



{t 



2\5 



G 



48 



logt' 



where the G-function of B2 is also of logrithmic type. However, in contrast with the 
Benney hierarchy, the second Poisson brackets of the dispersion expansion coincide 
with those of DZ brackets only up to e° 



15]. This inconsistency also pointed out in 



|20| by considering the commuting flows. 



VI. CONCLUSIONS 



We have studied several interesting properties associated with the bi-Hamiltonian 
structure of the Benney hierarchy. Starting with the Poisson brackets of hydrody- 
namic type we obtain the structure coefficients of an associative algebra character- 
izing the associated Frobenius manifold. This implies that there exist a function 
which generates the structure coefficients and this function in fact can be viewed 
as the genus zero primary free energy of a TFT. It turns out that the topological 
correlation functions at genus zero can be constructed by using the LG formula- 
tion. After appropriately define the two-point correlation functions as the Benney 



variables, the genus zero topological recursion relation of the TFT turns out to be 
the Benney equations. Furthermore, we use the twistor construction to derive the 
string equation associated with the TFT, which describe the dynamics of correlation 
functions in the full phase space. Finally, based on the approach of Dubrovin and 
Zhang we obtain the genus one correction of the Poisson brackets. We show that 
the same result can be reached by analyzing the OS brackets of the KB hierarchy 
in semi-classical limit. 

Inspite of the results obtained, there are some interesting issues deserve more 
investigations. 

• One knows that there exists a Legendre-type transformation between the free 



energy ( p.8|) and that of topological CP^ model |jT8|. In fact, one can check that 
the DZ brackets of the Benney hierarchy and those of dToda system in CP^ 
model [0,|3] coincide only up to e° via the Legendre-type transformation. 
A natural question is: Can we find a "quantized" version of the Legendre 
transformation between these two dispersionless hierarchies so that their DZ 
brackets match up to (? > 1? 

In the deviation of the string equation we borrow the method of twistor con- 



struction to obtain the result. In [1^, using recursion procedure, Dubrovin 



and Zhang establish the Virasoro constraints of genus zero for arbitrary Frobe- 
nius manifold. Then it would be interesting to know whether the solution of 
string equation of the Benney hierarchy satisfies Virasoro constraints using 
the methods developed in [^,^. Also, we wonder that whether there exists 
a matrix model associated with the Benney hierarchy such that the large-A^ 
limit of that reproduces the genus expansion of the hierarchy flows. 

Recently, the genus-two free energy T2 of the Benney hierarchy has been 
obtained in |T^ (see also fl^) by combining genus-two topological recursion 



relations and Virasoro constraints |T^. So we should extend the expansions 



( |5.8|) and (|5.9|) to O(e^) when the genus-two correction is included. Then it 
is quite interesting to see whether the expansions ( ^.lOj ) can be extended so 
that the OS brackets are matched with the DZ brackets up to O(e^) after 
appropriate differential substitutions. 

We hope to report on these issues in a forthcoming paper. 
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APPENDIX A: PROOF OF THEOREM 1 



Proof. Let us first derive the canonical Poisson relations 
we have 



By the chain rule 



dfif,,M) 


^fi^^M) \ 


( dp. 


dp 


dfi 


dM 


dp 
dM 


dX 


^g{^l,M) 


ag{^l,M) 


dM 


a/x 


dM I 


\ dp 


dX 



df(pM) 

dp 
dgip,M) 

djl 



dfip,M) 

dM 
dg{p,M) 

dM 



dp djl 

dp dX 

dM dM 

dp dX 



(Al) 



Taking the determinant of both hand sides of this equation and using relations ( |4.7] ) 
we get 



{/i,M} = {/i,M}. 



(A2) 



One can calculate the left hand side as 
d/j dM dM d/j 



{^^,M} 



dp dX dp dX 




dii , ^ , ^ dw,iT,f) _ 

'n(T,f) fixed ^"^ i=i '^^ 

dfi 

„(r,f) fixed '^P 
= 1 + (negative powers of p) 

where we have used the fact that the terms containing ( ^ j _ j^ , in the last 

*^ V 9M/«;„(r,T) fixed 

line cancel. Similar calculations can show that the right hand side contains only the 
non- negative powers of p. Therefore the both hands of ( |A^ ) are equal to 1, which 
proves the canonical relations ( [4.4|) . 

The Lax equations with respect to T„ are proved as follows. Differentiating 
equations (f4.8|) by T„ gives 



9f(p,M) 

dp 
dgjpM) 

dp 



dfjpM) 

dM 
dg{p,M) 

dM 



dT„ 



Using (|A1|), we can rewrite (^) as 

/dp. dp_ \-'^ / dp_ ^ 

\ dil dM_ ] dl3_ 

\ dp dX J \ dT„ , 



df{p,M) 

dp 
dg{p,M) 

dp 



dp djl 

dp dX 

dM dM 

dp dX 



df{p,M) 

dM 
dg{p;M) 

dM 



-1 





(A3) 



(A4) 



Since the the determinants of the 2x2 matrices on both sides are 1, the inverse 
can also be written explicitly. In components, thus, the above matrix ( |A4|) relation 
gives 



dM dfi dfi dM _ dM dfi dfi dM 

'dxWn ~ dXdf^ ~ 'dXdfn ~ dX dT^ 

dM (9/i (9/i dM _ dM djl djl dM 

'dpdf^ ~ dpdT^ ~ 'dpdT^ ~ dpdT^' 



(A5) 
(A6) 



The left hand sides of equations { \Kq ) and (|A6|) can be calculated just as we have 
done above for derivatives in {p,X). For the equation of ( |A5| ) , 

dM dfi dfi dM 



dX dTr, dX dTn 

dM\ ^/^ , 1 , ^dw,{T,f) _ 

/^ + ^ + -^ 7)Y — ^ 

i{T,f) fixed '^^ "-1 ^^ 



dfi 



i=l 



dfi 



dfi 
dX 



'dM' 



•^/^ J Wi{T,f) fixed ^-^' 



^^ +n;."- 



^dw,{T,f) _ 



and terms containing (^) - n ^ cancel. Thus, 
^ \dfj. Jwi{T,T) fixed ' 



dM dfi dfi dM difi""] 



>i 



dX dT„ dX dT„ 



dX 



+ (powers of p <0). 



(AT) 



Similar calculations show that the right hand side of ( [A5| ) contains only positive 
powers of p. Therefore only powers of p > 1 should survive. Hence 



dM dfi dfi dM d{ff' 



)>i 



dX dTn dX dTn 
For the equation (|A6|) , we have similarly 

dM dfi dfidM d{fi'')>o 



dX 



dB^ 
dX 



{M 



dp dTn dp dTn 



dp 

d(y%i 

dp 



+ (negative powers of p), 
+ (negative powers of p). 



Similarly, noticing the partial derivative ^, we can also show that the right hand 
of ( |A6| ) has Laurent expansion with only non-negative powers of p. Hence only 
nonnegative powers of p should survive. Thus 



dM dfi dfidM _ d{fi'')>i _ dB^ 
dp dTn dp dTn dp dp 



(A9) 



Using 



{fi,M} = {fi,M} 



the equations ([A8|) and ( |X9|) can be readily solved: 

d^ ^ _dji^dRn_ dji_dR^ ^ 

dTn~ dpdX ^ dX dp -^^"'^^' 

m _ _dM_dB^ mdB^ _ 

dTn dp dX ^ dX dp ^^-'^'>' 

which is nothing but T„-fiow part of the Lax equations ( [4.2| ) and ( [4 .61 ) . The T„-fiow 
part of the Lax equations can be proved in the similar way. This completes the 
proof of the theorem. 



APPENDIX B: A DEVIATION OF G-FUCNTION 



For the quasi- homogeneous primary free energy ( p.8| ), one can introduce the 
canonical coordinates u^,v? determined by 

dei{g'^ it) - ur]'^ it)) = 

which gives 

u^ = e - 2Vt2, u^ = t^ + 2Vt2. 

Then the topological metric 77*-' and its inverse rjij in the canonical coordinates are 
given by 



It has been shown []T5[ that the G-function can be expressed in the following formula 



G{t\t') = \og^ 



24 



where J is the Jacobian of the transformation from the canonical coordinates to the 
flat ones, which is easily obtained as 

and Tj is the tau-function of a solution in the theory of isomonodromy deformation 
defined by the quadrature |]26[ 



d\ogTi = Y,Hidu' (B2) 

where the quadratic Hamiltonian Hi has the form 



1^ V^- 






with 



Vij = -{u -u-')'jij{u), 'jij{u) = 



VjjW 



_d_ 

du'- 

From flBll) we have 



Note that di = -^ and Vij = —Vji. 



i.Au) = 1 !^ !;!V!:i ) , v., = {i 7 



M^-U'^) 2(«2_„1) \ ,, / 

-i 1 

2(«2_„i) 2(«2_„i) / V 2 



which together with ([B2|) imphes 

n = in' - uY'^', 
and hence, up to an additive constant 

^,2 _., 1^-1/8 1 
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